Energy shifted level densities in the rare earth region 
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The density of accessible levels at low spin in the ('^He,Q:7) reaction has been extracted for 
the ^^^'^^^Dy and I'^^'^^^Yb nuclei. The energy shift between the level densities of the even-odd 
and even-even isotopes is measured as a function of excitation energy. The results are compared 
with predictions from various semi-empirical models. The energy shift procedure works well for 
excitation energies between 3.5 and 7 MeV in the even-even nucleus, provided that a proper level 
density function is used. The experimental energy shift is close to the pairing gap parameter A. 
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I. INTRODUCTION 



The Fermi gas model for finite nuclei has been shown 
to account for the nuclear level density at high excitation 
energy |jl|. The model describes the nucleus as a gas of 
non-interacting fermions confined to the nuclear volume 
and neglects shell effects and pair correlations. 

The level density of odd-odd nuclei is relatively high 
already around the ground state due to the coupling of 
the two valence particles to the core. However, for odd- 
mass and even-even nuclei the level density is consid- 
erably lower due to the formation of Cooper pairs. A 
handful semi-empirical approaches has been suggested to 
describe this effect of pair correlations by a simple energy 
shift of the level density function. 

In the conventional shifted Fermi gas model |p|-0| the 
excitation energy is shifted using the pairing energy pa- 
rameter A. The shifts are A and 2 A for odd-mass and 
even-even nuclei, respectively, yielding approximately the 
level density found for the neighboring odd-odd system. 
This description turned out to be too rigid to reproduce 
the level densities at low and high excitation energies, 
simultaneously. A two-component level density formula 
with energy shifts was later introduced [Hj|]. Here, the 
first '^ 10 MeV of excitation energy is described by a 
constant temperature formula, and at higher energies the 
shifted Fermi gas model is applied. 

A simpler and rather well working version is the back- 
shifted "'^ Fermi gas model H , where the Fermi gas formula 
is used for all excitation energies. The model has only 
two parameters: the back-shifted energy and the level 
density parameter a, both being free parameters in order 
to fit the data. 

There are several unclear points in using these ap- 
proaches. The main questions concern the functional 



form of the level density and the justification of a shift 
of the excitation energy to describe the level densities of 
neighboring nuclei. The fact that the extracted shifts sel- 
dom coincide with the paring gap parameter A (or 2A) 
indicates that one or both assumptions are not fulfilled. 

In the vicinity of the ground band, levels can be 
counted reliably up to a certain excitation energy, typi- 
cally 1.5 MeV in even rare earth nuclei. The level density 
can also be derived by the level spacing of neutron res- 
onances at the neutron binding energy _B„. In between 
these energies few experimental results are available. In 
the A = 40 - 60 mass region, fiuctuation widths and 
charged particle emission spectra have been employed. 
Recently g], the level density of ^^As and ™Ge was mea- 
sured using the (^^C,p) compound reaction. The data 
cover 5-24 MeV of excitation energy, but eventual fine 
structures are smeared out due to the high average spin 
and broad spin distributions in such studies. 

The Oslo group has extracted level densities from mea- 
sured 7-ray spectra p|-p^. The method, which is de- 
scribed in Ref. ||l^, allows for the simultaneous extrac- 
tion of the level density and the 7-strength function over 
a wide energy region. In these studies, the transferred 
spin in the (■^He,Q;) reaction is approximately 2 - 6 ?i, 
and the nuclear system is believed to thermalize prior to 
7-emission. 

The subject of this work is to extract experimental en- 
ergy shifts using the method of Ref. |T^ and to investi- 
gate the quality of the energy shift procedure as function 
of excitation energy. Furthermore, it is interesting to 
compare the value of the energy shift parameter to the 
pairing gap parameter A. 



* E-mail: magne.guttormsen@fys.uio.no 

'^For light nuclei the ground state position has to be shifted 
to lower energies compared to the fit performed with the 
conventional shifted Fermi gas model; therefore the notation 
"back-shift". 



II. EXPERIMENTAL METHOD 

The experiments were carried out with 45 MeV '^He- 
projectiles at the Oslo Cyclotron Laboratory (OCL). 
The particle-7 coincidences are measured with the CAC- 
TUS multidetector array [Q using the {'^lle,a^) reaction 



162,163 



Dy and 



172,173 



Yb self-supporting targets. The 



on 

charged ejectiles were detected with eight particle tele- 
scopes placed at an angle of 45° relative to the beam 
direction. An array of 28 Nal 7-ray detectors with a to- 
tal efficiency of ~15% surrounded the target and particle 
detectors. 

The experimental extraction procedure and the as- 
sumptions made are described in Refs. P,13[ and ref- 
erences therein. The level density is deduced from 7-ray 
spectra recorded at a number of initial excitation ener- 
gies E, determined by the measured a-energy. These 
data are the basis for making the first generation (or pri- 
mary) 7-ray matrix, which is factorized according to the 
Brink- Axel hypothesis ||l^,|l^ as 



PiE, E^) ex p{E - E^)(7{E^), 



(1) 



where the level density p and the 7-energy dependent 
function a are determined by an iterative procedure. 

To obtain the level density, the first trial function for p 
is simply taken as a straight line and the corresponding 
a is determined by Eq. (1). Then a x^ minimum is calcu- 
lated for each parameter, keeping the others fixed. This 
procedure is repeated about 50 times, until a global least 
square fit to the ~ 1400 data points of P is achieved. 

It has been shown p3| that if one solution for p and a 
is found, functions of the form 



and 



p{E - E^) -^ Aeyi^[a{E - E^)]p{E - E^) (2) 



(j{E^)~^ BeyiY>{aE^)a{E^), (3) 



give exactly the same fit to the P{E,E.y) matrix. The 
values of A, B and a can be determined by additional 
conditions. The A and a parameters are used for abso- 
lute normalization of the level density p: They are ad- 
justed to reproduce (i) the number of levels observed in 
the vicinity of the ground state and (ii) the neutron reso- 
nance spacing at the neutron binding energy i?„. Further 
details on the method, the iteration procedure and the 
simulation of errors are given in Ref . ||l^ . 

In the following we will not discuss the 7-energy de- 
pendent function cr, but concentrate on the level density 
P- 



III. EXPERIMENTAL AND SEMI-EMPIRICAL 
LEVEL DENSITIES 

The experimental level densities for the ^^^•^^'^T)y and 
i7i,i72y]3 nuclei are shown as data points in Figs. 1 and 
2. In the extraction technique, we exclude data with 7- 
energies below 1 MeV due to methodical problems in the 
first generation spectra. Therefore, the level density is 
generally determined only up to E = Bn — 1 MeV. Re- 
cently ||l^ , thermodynamical aspects have been discussed 
from the very same level density curves as shown here. 
The figures include the level densities (solid drawn lines) 
obtained from counting known levels [18]. These known 
densities are seen to agree with our experimental data 
points. The level densities from counting are reliable up 
to an excitation energy of E ^ 0.8 MeV and 1.7 MeV for 
the odd- and even-mass cases, respectively. Above these 
energies, more than 100 levels are present per MeV and 
it is experimentally difficult to resolve all levels. There- 
fore, the density of known discrete levels drops at higher 
excitation energies. 

The level densities for the ^^^Dy and ^^^Yb isotopes are 
about five times higher than for the neighboring ^^^Dy 
and ^^^Yb isotopes. The latter isotopes seem to exhibit 
the same slope at high excitation energy. However, the 
presence of bumps modifies this simple picture, in par- 
ticular at low excitation energies. 

The energy region up to '^ 5 — 10 MeV has been de- 
scribed by the constant temperature formula [5,6] given 

by 



p = Ccxp{E/T), 



(4) 



where the normalization factor C and the temperature T 
are constants. Also level densities based on the Fermigas 
model are frequently adopted in this energy region [^-Q] 



P 



exp[2> 



12V2ai/4[/5/4^' 



(5) 



where a is the spin cut-off parameter and U is the 
back-shifted energy. As examples of such approaches 
Figs. 1 and 2 also include level densities from Gilbert 
and Cameron [5] (dashed curves) and from von Egidy et 
al. [6] (dash-dotted curves). Full details on the formulae, 
parameterizations and choice of parameters are given in 
Refs. [5,6]. 

The level densities of Gilbert and Cameron are de- 
scribed by Eq. (4) in the excitation region below ~ 5 
MeV, and at higher energies they use Eq. (5). The de- 
scription (dashed curves) is rather poor, except for ^^^ Yb. 
Gilbert and Cameron give temperatures that are lower in 



^The temperature T can be extracted from the inverse of 
the slope of the logarithm of the data points, i.e., 1/T = 
dlnp/dE. 



the even-even systems, contrary to the tendency of our 
data^. It seems that these authors have, for the even- 
even nuclei, anchored their constant temperature level 
density curves to the ground state band, rather than to 
levels at ~ 1.7 MeV, where the two quasiparticle regime 
appears. Also scarce data at the time their compilation 
was made (1965) could be the reason for the poor agree- 
ment. Even so, we think the two-component level density 
is a reasonable approach. In the first MeV of excitation 
energy, nucleon pairs (Cooper pairs) are broken and thus 
prevent the temperature to rise as fast as predicted by 
the Fermi gas formula. This mechanism is discussed in 
Refs. [11,12], and references therein. For excitation en- 
ergies around and above the neutron binding energy, the 
Fermi gas conditions are probably fulfilled. Here, the 
pairing correlations are quenched and a high density of 
single particle levels is present. 

von Egidy et al. [6] have tested both the constant tem- 
perature and the back-shifted Fermi gas formulae in this 
region. They find that both approaches give similar ^- 
fits to experimental data. The suggested temperatures 
are close to 0.6 MeV for all four nuclei, almost 0.1 MeV 
higher than our data indicate. In Figs. 1 and 2 we show 
only the Fermi gas results (dash-dotted curves). Within 
this model both the level density parameter a and the 
correction to the back-shift C\ {U = E — A — Ci) are 
based on global paramctcrizations as function of the mass 
number A. With this restriction, one may say that the 
level density curves describe our data points rather well. 
However, a clear shortcoming is that these expressions 
increase too slowly as function of excitation energy. The 
experimental level densities, in particular for the Yb iso- 
topes, show a functional form closer to the constant tem- 
perature formula"^. 

Global fits for all mass numbers can give deviations 
of up to a factor of 10 from known average neutron res- 
onance spacings [6]. Of course, better local fits to ex- 
perimental data could be achieved, both with Eq. (4) 
and/or Eq. (5). However, a common approach is difficult 
to construct since all four nuclei exhibit different func- 
tional forms. These variations are probably connected 
with details in the quenching of the pair correlations in 
the individual nuclei. 

In the following, we will not propose any new semi- 
empirical formulae, but rather focus on the feature that 
the level densities for the neighboring isotopes exhibit 
similar shapes as function of excitation energy. 



IV. COMPARISON BETWEEN EXPERIMENTAL 
AND SEMI-EMPIRICAL ENERGY SHIFTS 

It is commonly believed that neighboring odd-odd, 
odd-even, even-odd and even-even isotopes reveal the 
same level density if a proper shift is applied to the ex- 
citation energy. This is indicated in Fig. 3, where the 
energies are compared to the ground state of the odd- 
odd nucleus. The even-even nucleus has its ground state 
pushed down by Ap -I- A„ relative to the odd-odd system 
due to pairing interactions. The odd-mass isotopes rep- 
resent cases in between these limits. With the present 
experimental data, we have the opportunity to test how 
well the energy-shift procedure works for the level densi- 
ties of even-odd (eo) and even-even (ee) systems. 

Neglecting collective excitations and two-body forces, 
various models can describe the level density of the odd- 
odd (oo) nucleus rather successfully, e.g. the Fermi gas 
model. The level density of the other neighboring nuclei 
can then be estimated by 



PociE) =^ PooiE ~ A„), 

Pco{E) = Poo{E - Ap), 

Pcc{E) = Poo{E ~ A„ - Ap), 



(6) 

(7) 
(8) 



where E is the excitation energy. 

The pairing gap parameters Ap and A„ can be deter- 
mined from empirical masses of a sequence of isotones or 
isotopes where [Q 

Ap = ^\Sp{N, Z + 1)- 2Sp{N, Z) + Sp{N, Z - l)\, (9) 

A„ = i|5„(7V + l,Z)- 2Sn{N, Z) + Sn{N - 1, Z)\, (10) 

and iSp and Sn are proton and neutron separation ener- 
gies |18| , respectively. The pairing gap parameter can 
alternatively be calculated by the empirical formula ]1^ 

A = l2A-^/'^MeY, (11) 

which is valid for both neutrons and protons. 

Equations (9) and (10) depend on the proton (Z) and 
neutron [N) numbers and should in principle give the 
best estimate. However, Eq. (11) gives a smooth function 
which neglects local shell effects, and is probably more 
correct if A is interpreted as a pure pairing parameter. 

From the extracted level densities for the I6i,i62j-)y ^^-^^^ 
i7i,i72y|3 nuclei, we can investigate the energy shift nec- 
essary to apply in order to simulate the level density in 
neighboring even-odd and even-even isotopes. The en- 
ergy shift 5{E) is defined as the necessary shift of the 
even-odd nucleus level density in order to describe the 
level density in the neighboring even-even nucleus 



^Gilbert and Cameron obtain a poor constant temperature 
fit for ^^^Yb, probably because they try to fit the level density 
of the ground band. 



pUE)^P.o{E-5{E)). 



(12) 



V. CONCLUSION 



In Fig. 4 the resulting 5{E) curves are plotted as func- 
tion of the excitation energy E measured in the even-even 
nucleus. In the excitation energy region between 3.5 and 
7 MeV the energy shift is rather constant giving 5 = 
1.13(7) and 0.84(10) MeV for I6i,i62j3y ^^^^ i7i,i72Yb^ 
respectively. 

The corresponding (5-values should also compare with 
the effective pairing gap parameters of 

A^ff(i62Dy) = Ap(i62Dy) + A„(i62Dy) - ApC^^By) 

(13) 



and 



^eff (i72Yb) ^ ApC^^Yh) + A„("2Yb) - Ap(i^iYb), 



(14) 



where we apply Eqs. (9) and (10) for the pairing gap pa- 
rameters. In Table 1 these values (and the pairing gap 
parameters calculated from Eq. (11)) are compared to the 
experimental 5-values. These values coincide rather well 
within less than 0.2 MeV, and Fig. 4 shows that A'''^ is in 
good agreement with the observed energy shifts in the 3.5 
- 7 MeV excitation region. The agreement is not that im- 
pressive when comparing experiment with energy shifts 
obtained from Eq. (11) and with semi-empirical level den- 
sity functions. The shifts from Gilbert and Cameron [5] 
(dashed curves in Fig. 4) deviate strongly from the ex- 
perimental data, as also indicated from the ^-values'' of 
Table 1. The shifts from von Egidy et al. [6] (dash-dotted 
curves) are determined by the experimental pairing gap 
A and the slow varying back-shift correction Ci. There- 
fore, these shifts coincide almost excactly with the A°*^ 
values. The small deviations seen in Fig. 4 are due to 
the 0.1 MeV""'^ increase in the level density parameters 
for the even-even systems. 

Both the two-component formula of Gilbert and 
Cameron [pi as well as the Fermi gas formula of Ref. R] 
give 5- values deviating with about 0.3 - 0.5 MeV (see Ta- 
ble 1). This is probably due to the free adjustment of S 
and other parameters, and indeed the shifts have been as- 
sociated with large uncertainties by these authors. The 
role of S in this type of approaches is not a pure en- 
ergy shift, but may also include a compensation for the 
unphysical form of the adopted analytical level density 
function. The same conclusion is evident from the com- 
pilations of Refs. |||-|^, where the extracted energy shifts 
scatter typically within ± 0.5 MeV in this mass region. 



The shifting of excitation energy in order to simulate 
the level density of neighboring isotopes works well us- 
ing realistic level density functions. The level densities 
follow each other rather close as function of excitation 
energy in the 3.5 - 7 MeV region, and the energy shifts 
coincide within 0.2 MeV with the pairing gap parameter 
A. Of the approaches studied here, only the formalism 
of von Egidy et al. has this feature included. Below 3.5 
MeV of excitation energy, nuclear structures assigned the 
various nuclei prevent the use of a simple energy shift 
procedure. In particular, the even-even isotopes reveal 
bumps in the level density function due to the breaking 
of Cooper pairs. 

The application of analytical level density functions 
may give energy shifts that deviate significantly from the 
experimental pairing gap parameters. In this work we 
have demonstrated that the functional form and/or the 
parameters used are not fully appropriate for the descrip- 
tion of level densities below 7 MeV of excitation energy. 

Probably, no simple level density formula can describe 
simultaneously the four nuclei investigated here. The 
Yb isotopes exhibit a constant temperature-like behav- 
ior, while the Dy isotopes are closer to the back-shifted 
Fermi gas prediction. Nevertheless, we find that the pa- 
rameters used for the semi-empirical formulae should un- 
dergo a revision. Here, all new low-lying levels should be 
included together with recent information on resonance 
level spacings. This effort, combined with a refined two- 
component formula, like the one of Gilbert and Cameron, 
could probably give better analytical formulae for future 
use. The analytical expressions should have the ability 
to give a constant energy shift between the level densities 
of neighboring isotopes, as observed in this work for the 

161,162i3y ^^^ 171,172Yb nuclei. 
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Table 1: Energy shift 5 extracted between the even-odd and even-even isotopes". 



Parameter (5 or A (MeV) 


i«2Dy 


172Yb 


5 from present data 


1.13(7) 


0.84(10) 


A"^** from separation energies, Eqs. (13,14) 


1.05 


0.93 


A from empirical formula, Eq. (11) 


0.94 


0.91 


5 from back-shifted Fermi gas [7] 


0.88(50) 


1.15(50) 


& from two-component level density [5]: 






Energies below ~ 5 MeV, Eq. (4) 


0.81(30)'' 


1.30(30) 


Energies above ~ 5 MeV, Eq. (5) 


0.70(20) 


0.69(20) 



° The energy shifts from semi-empirical level densitiy formulae are constants based on parameters used in Rcfs. 
The energy shifts of von Egidy et al. [6] are not listed since they coincide with A"**^. 
^ The shift is calculated from the i^O-i^ij^y parameter sets. 



[5-7]. 
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Excitation energy [MeV] 



FIG. 1. Observed level densities for '^ Dy as functions of excitation energy. The experimental data points are compared 
to the density of known levels at low excitation energy (solid lines) . The figure also includes the semi-empirical level density for- 
mulae of Gilbert and Cameron [5] (dashed curves) and von Egidy et al. [6] (dash-dotted curves). Upper and lower points/curves 
are for ^^^Dy and ^^^Dy, respectively. Since Gilbert and Cameron give no parameters for ^^^Dy, we use the ^^°Dy parameter 
set. 
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Excitation energy [MeV] 



FIG. 2. Observed level densities for '^ Yb as functions of excitation energy. The experimental data points are compared 
to the density of known levels at low excitation energy (solid lines) and the semi-empirical level density formulae of Gilbert 
and Cameron [5] (dashed curves) and von Egidy et al. [6] (dash-dotted curves). Upper and lower points/curves are for ^^^Yb 
and ^^^Yb, respectively. 
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FIG. 3. Illustration of how the level density for various nuclei (oo, oe, eo, and ee) can be estimated by proper energy shifts. 
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FIG. 4. The observed energy shift S between Dy and ^ Dy (upper part) and Yb and ^ Yb (lower part). The pairing 
gap parameters A*' , evaluated from Eqs. (13) and (14), are displayed as solid lines for comparison. The energy shifts obtained 
from Gilbert and Cameron [5] (dashed curves) and von Egidy et al. [6] (dash-dotted curves) are also displayed. 



